YGHP-03-31 



Improved gauge action on an anisotropic 

lattice II 

- Anisotropy parameter in the medium 
coupling region - 

S. Sakai a and A. Nakamura b 



a 



Faculty of Education, Yamagata University, Yamagata 990-8560, Japan 
h RIISE, Hiroshima University, Higashi- Hiroshima 739-8521, Japan 



Abstract 

The quantum correction of the anisotropy parameter, n, is calculated for £ = 2 
and 3 in the (5 region where numerical simulations such as hadron spectroscopy are 
currently carried out, for the improved actions composed of plaquette and rectan- 
gular 6-link loops. The (5 dependences of rj for the renormalization group improved 
actions are quite different from those of the standard and Symanzik actions. In 
Iwasaki and DBW2 actions, n stays almost constant in a wide range of (5, which 
also differs from the one-loop perturbative result, while in the case of Symanzik 
action, it increases as (3 decreases, which is qualitatively similar to the perturbative 
result, but the slope is steeper. In the calculation of the rj parameter close to and 
in the confined phase, we have applied the link integration method to suppress the 
fluctuation of the gauge fields. Some technical details are summarized. 



1 Introduction 



Anisotropic lattices, with temporal lattice spacing smaller than the spatial 
one, provide an effective method of precise Monte Carlo calculations of, for 
example, heavy quark systems, glueball masses, finite temperature properties 
of the QCD. Properties of the anisotropic lattice with the standard plaquette 
action have been studied by several groups [1,2, 3]. 

On the other hand, improved actions are proposed to obtain numerical 
results close to the continuum limit on relatively coarse lattices. Therefore it 
is worth studying the anisotropic properties of the improved actions. 
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In the previous paper, we studied the anisotropic lattice for a class of im- 
proved actions in the weak coupling region, mainly using the perturbative 
method[4]. The improved actions we considered are written in terms of the 
plaquette and 6-link rectangular loops as 

S oc E [C P(1 x 1)^ + dP(l x 2U , (1) 

where C and C\ satisfy the relation C + 8C1 = 1. The improved actions 
frequently used in the simulations correspond to the following parameters: 
Ci = —1/8 (Symanzik's improved action[6]), C\ = —0.331 (Iwasaki's improved 
action[5]), and d = -1.4088 (QCDTARO's DBW2 action[7]). 

For these types of actions, we can formulate the anisotropic lattice in the 
same way as the standard plaquette action, 

s a = &(-!- E E Pa + E E p«), (2) 

SB a; j>j x 



where = ^ f3 a f3 T1 and is a bare anisotropic parameter which controls the 
anisotropy in the space and time directions. The anisotropy is defined as the 
ratio of the lattice spacing in the spatial (a a ) and the temporal (a T ) directions, 
0? = a a /a T . 

Due to quantum correction, ^ R is not equal to £g. Therefore it is essential 
to know their relationship before large scale simulations on the anisotropic 
lattice with improved actions are carried out. 

The effects of quantum correction have been studied using the parameter rj 
defined by 

n = I 5 - (3) 

In the weak coupling region, it is calculated perturbatively[l],[4]. The one- 
loop perturbative results have been very impressive in the sense that as —C 1 
increases, qualitative change is observed in the behavior of i] as a function of 
(5. If it is parametrized as 

iV 

i7(e,0,Ci) = 1 + jifctf, CO, (4) 

the coefficient 771 decreases as —C\ increases, and at around C\ ~ —0.18, it 
reaches to zero and then becomes negative. Therefore r\{0) of Iwasaki action 
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and DBD2 action decreases as (3 decreases, while in the case of the standard 
and Symanzik actions, it increases. Namely, they have opposite (3 dependences. 



The natural question is what would be the behavior of 77 in the smaller (3 
region, where the perturbative calculation breaks down. In this work, we will 
focus on the rj parameter at £r = 2 and 3, and calculate it in the intermediate 
(3 region where most current numerical simulations are carried out. 



In Section 2, we discuss appropriate regions of (3 for the improved actions to 
evaluate 77, and explain some details of the calculation: matching of the lattice 
potential in the spatial and temporal directions, and the method of eliminating 
the effects of self-energy terms of lattice potential on an anisotropic lattice. For 
the standard action, the effect of the self-energy terms on 77 has been studied 
by the Bielefeld group [3] . They reported that the effects on 77 are small about 
1% throughout the parameter range that they studied. We discuss the effects 
for improved actions and show, in Section 3, that they are not large. 



In Section 3, 77 for each improved action is presented. The 77 behavior in 
the intermediate (3 region is quite different for each improved action. Iwasaki 
and DBW2 actions are qualitatively different in the one-loop perturbative 
calculation. 



Section 4 is devoted to the discussion and conclusion. It is found that for 
Iwasaki action, 77 is close to unity in the region (3 > 2.5, indicating that 
the detailed calibration of for a given £ R is not important unless a very 
precise simulation is carried out. For DBW2 action, 77 remains essentially con- 
stant, indicating that rough calibration gives a reasonable estimation of 77. For 
Symanzik action some detailed calculation of 77 is necessary, as in the case of 
standard action. 



For the calculation of the 77 parameter, measurements of the large Wilson 
loops are required. Large Wilson loops suffer from huge fluctuation of the gauge 
fields, particularly in the confined phase. To suppress the fluctuation, the link 
integration method has been proposed [8,9,10]. In this work, we have applied 
the link integration method in small (3 regions. Here, it is very important to 
choose an adequate radius (optimal radius) of integration in the complex plane. 
For the standard action, the optimal radius was studied by the Bielefeld Group 
[3]. In appendix A, we describe the optimal radii for Symanzik and Iwasaki 
actions. 
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2 Calculation of the anisotropy quantum correction, rj 



2. 1 Region of coupling constant to be studied 



In this work, we calculate the rj parameter in the regions of (3 where most 
numerical calculations are currently carried out. In the case of standard action, 
hadron spectroscopy in the quench approximation has been reported for 5.7 < 
j3 < 6.2 [11]. In these coupling constant regions, the light hadron masses 
are reproduced within 10% accuracy, which may be the limit of the quench 
approximation. Therefore, we calculate the n parameter around these lattice 
spacings for the improved actions. 

In order to estimate the lattice spacing for the improved actions, we use the 
critical f3 (flcrit) of finite temperature transition. For standard action, (3 = 6.05 
corresponds to the finite temperature transition point of the N T = 8 lattice [12]. 
Therefore we estimate Pent at N T = 8 for the improved actions. For tree level 
improved Symanzik action, (3crit s are reported for = 3, 4, 5 and 6[13], and 
for Iwasaki action, they have been calculated at N? = 4, 6 by the Tsukuba 
group[14] and at Nt = 8, by the Yamagata-Hiroshima collaboration [15]. For 
DBW2 action, they are reported by the QCDTARO collaboration [7] for N T = 
3, 4 and 6. 

Pent at Nt = 8 is estimated using the two loop asymptotic scaling relation 
for lattice spacing, 

a = A { T ] exp( "W' (5) 



where a is the lattice spacing and c = 6i/(26q), with bo = ll/(47r) 2 and b\ = 
102/(47r) 4 . We applied two methods for the estimation of Pent at Nt = 8. In 
method 1, we use (3 C rit at N T = 6 of the same action to estimate (3crit(N T = 8) 
using T c = l/(aN T ). In method 2, we use the Pcrit(N T = 8) of the standard ac- 
tion to evaluate ficritiNr = 8) of the other action with the A parameter [4, 16]. 
They are summarized in the table 1. 

In the case of Symanzik action, the estimations of ficrit by the two methods 
Table 1 

Estimation of ficrit at Nt = 8 for various actions 





Pcrit{data) 


Method 1 


Method 2 


Minimum (3 


Standard 


6.05(iV T = 8) 








Symanzik 


4.31(iV r = 6) 


4.57 


4.56 


4.5 


Iwasaki 


2.52(iV T = 6) 


2.78 


2.32 


2.5 


DBW2 


0.936(iVr = 6) 


1.28 




1.1 



coincide with each other. For Iwasaki action, some discrepancy is observed 
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between the two estimations. The first method gives a closer result to that of 
Ref.[15], in which f3 Crit = 2.73 ~ 2.75. 

For DBW2 action, (3 C rit estimated using the A ratio becomes negative. In 
this (3 region, the deviation from the perturbative scaling relation is quite large 
for this action. Therefore, for the estimation of f3c r ti(^T = 8), we plot Pent 
at N t = 3,4 and 6 and simply extrapolate it. It becomes about (3 = 1.1 with 
large ambiguity. Therefore we will study the action until (3 = 1.0. 

In the table 1, we show also the minimum of f3, for which we calculate the 
rj parameters. 

2.2 Subtraction of the self-energy contribution from the lattice potential 

The renormalized anisotropy £r is defined as the ratio of the lattice spacings 
in the spatial and temporal directions, £r = a a /a T . In the quench approxima- 
tion, the lattice potential has been used as the probe of lattice spacing, and 
is defined as the ratio of Wilson loops, 

The lattice potentials in the spatial (V s ) and temporal directions (Vj) are deter- 
mined by the Wilson loops in the space-space plane and the space-time plane, 
respectively. The lattice potential defined in Eq.6 will become independent of 
the position p when p becomes large where the lattice artifact disappears. In 
the following, we will discard p in the artifact-free region, but write it explicitly 
when we discuss its effect. 

The matching of the potential in the spatial and temporal directions [2,17] 
has been used for the determination of £r as 

V a (tB,r)=V t (Z B ,t = Z R xr). (7) 

In the calculation of rj, we fix the renormalized anisotropy £r, and then search 
for the point of where the Eq.7 is satisfied [2]. Using these and £r values, 
the 7] parameter is calculated. 

The lattice potential defined by Eq. 6 suffers from the self-energy term. In 
this article, we assume the simplest parameterization for it: 

V s ^B,r) = V s °^ B ) + V s L ^ B ,r), (8) 
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where V S L is the lattice potential free of self-energy contributions. The temporal 
potential, V t (£e,t), is treated similarly. On an anisotropic lattice, V^^b) an d 
K°(£b) ma y be different from each other due to the anisotropy. Therefore the 
matching of the potential should be applied for the self-energy-free parts, V S L 
and V t L . 

In order to eliminate the effect of the self-energy term V°, we define the 
subtracted potential: 

V s Sub (Z B ,r,r ) = V s (£ B ,r) - V a (t B ,r ) = V s L (£ B ,r) - V s \i B ,r ). (9) 



V t Sub is defined in a similar manner. 

The subtraction points ro and to are chosen to satisfy to — £,R r o, and the 
matching of the potential, V t L (t = ^rTq) = V s L (r ), should also be satisfied 
there. Namely, at r , the lattice potential should be free of the lattice artifacts. 
This condition is satisfied if r is large. At large r , however, the fluctuation of 
the potential increases, and simulations with high statistics on a larger lattice 
are required. Therefore, tq should be chosen to be as small as possible, where 
the lattice artifacts will be sufficiently small. 



2.3 Matching method 



As an example, we will show details of the determination of 77 in the case 
of Iwasaki action at f3 — 4.5, £r = 2 on the 12 3 x 24 lattice. 

Let us start with the determination of the subtraction point r . In order to 
reduce statistical error, small r is preferable. In the small r region, however, 
systematic error due to lattice artifacts becomes large. The optimal choice of 
r requires careful study by trial and error. 

First we calculate the ratio 

mB - p ' r) = m^J^7y (10) 



where V s and V t mean the lattice potential in the space and time directions, re- 
spectively, and include the self-energy contributions. Our results are displayed 
in the Fig. 1. As r increases, the ratio R(p, r) approaches an asymptotic value. 
It is seen that at r = 1, the deviation from the asymptotic value is rather 
large, which may be due to lattice artifacts. Therefore, we first choose r = 3, 
and calculate the subtracted potentials of Eq.9, and then use them in Eq.10 
to obtain ratio R. 
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1.04 



1.02 



R(p,r) at 5 =2.1 




1.00 



0.98 



2 3 r 4 



Fig. 1. Ratios given by Eq.10. 
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Fig. 2. £b dependences of the ratios for the subtracted potentials. 

The results for £ B = 2.0, 2.1 and 2.2 are shown in Fig. 2. The ratios are 
shown individually for each p and r. We proceed to look for the point where 
the ratios satisfy the relation R(p, r, £b) — 1- We fit the three points by the 
second-order polynomial of and find the solution 

R(p,r,tB)=c + c 1 ZB + c 2 g } = l. (11) 



The coefficients Co, c\ and c 2 are determined by the three data points of 



7 



2.2 



T|(P,0 



2.1 ■ 
2.0 



1.9 



1.8 



10 15 20 pxr 25 30 



Fig. 3. r] = £r/£b at each p and r, where £r (£#) stands for the renormalized (bare) 
anisotropy. 



r, £b). Using the solution of Eq.ll, the ratios i] = £r/£b are determined 
for each p and r and are shown in Fig. 3. 



In order to avoid the lattice artifact, we choose the data with p > 3 and 
r > 4. In this region, r)(p,r) is almost independent of p and r. rj at this (5 
is determined as the average of the data. The errors are estimated by the 
jackknife method; the data after thermalization is grouped into 10 blocks 
and they are used as independent data. In this way the result becomes rj = 
0.9766 ± 0.0044 when r = 3. 



The same analysis is carried out by taking r = 2. The result becomes 
rj = 0.9764 ± 0.0039. The result changes little in this case. However, if we 
choose r = 1, the results differ significantly from those of r =2 and 3. 
Analyses are carried out at other values of (3 and £. There are several cases in 
which a slight difference is observed between r$ = 2 and r$ = 3. Therefore in 
the case of Iwasaki action, we choose r = 3 for all values of f3 and £b in this 
work. 



We carried out the same studies for Symanzik and DBW2 actions. In these 
cases, the subtraction point becomes r = 4. This indicates that lattice arti- 
facts are larger for these actions. 
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3 Result for the quantum correction of the anisotropy, 77 

3.1 Simulation parameters, numerical results and self-energy contribution 

The simulations are mainly carried out on the 12 3 x 12£r lattice. For some 
values of (3 and £, a larger 16 3 x 16£ R lattice is employed in order to study the 
size dependence. It is found that the lattice size effect is small. 

The gauge configurations are generated by the heat bath method with over- 
relaxation [18,19]. Typical numbers of the Monte Carlo (MC) data for the 
calculation of p,r) are a few tens of thousand after thermalization of 

about 10 4 MC sweeps. However, as j3 decreases and approaches the finite 
transition point or goes into the confined phase, both the necessary number 
of MC data and the number of thermalization sweeps increase. For the calcu- 
lation of t?(£ — 2) at (5 — 2.5 for Iwasaki action, we used 1.5 x 10 6 data after 
thermalization of 3.5 x 10 5 MC sweeps. 

In order to suppress the fluctuation of the gauge fields in the calculation of 
large Wilson loops, we applied the link integration method [9,10]. It is used 
for the calculation of the lattice potentials at f3 — 2.5 and 2.56 for Iwasaki 
action and (5 — 4.5 for Symanzik action. Technical details will be presented in 
the appendix. Here we only notice that, in the case of improved actions, the 
effect of the link integration is reduced due to the presence of the rectangular 
6-link loops. 

Our results of 77 are summarized in the tables 2 to 7. In order to show 
the effects of self energy terms in the lattice potential, we have presented the 
results for 77, which are obtained without subtracting the self energy terms in 
the rj Nosub column in these tables. It is found that the difference between them 
is less than ~ 1% for Symanzik and Iwasaki actions. This is consistent with 
the result for the standard action obtained by Bielefeld group [3]. However, 
for DBW2 action, the difference increases. It amounts to a few percent but is 
still small. Therefore, except for the case of the simulation with a few percent 
accuracy, it is safe to use rj Nosub , which has been reported at the XVIIth 
International Symposium on Lattice Field Theory at Pisa (Lattice '99) [21]. 



3.2 Symanzik action 

In Fig. 4, we show the values of the 77 parameter as a function of (5 for 
Symanzik action. The qualitative behavior of 77 as a function of (5 is the same 
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Table 2 

rj = £r/£b for Symanzik action at £r = 2, where £r and £b are the renormalized 
and bare anisotropics, respectively. At (3 = 8.0 and 4.5 the simulation is carried out 
on 16 3 x 32 lattice to study the size dependences. They are shown in the line with 
the symbol *. 

Symanzik Action at £r = 2 



Table 3 

rj as a function of 



(3 


V 


N oSub 


10.0 


1.0227 ± 0.0097 


1.0271 ± 0.0031 


8.0 


1.0393 ± 0.0191 


1.0391 ± 0.0020 


6.0 


1.0381 ± 0.0097 


1.0500 ± 0.0029 


4.5 


1.0980 ± 0.0255 


1.1011 ± 0.0021 


8.0* 


1.0284 ± 0.0021 


1.0232 ± 0.0039 


4.5* 


1.1095 ± 0.0122 


1.1040 ± 0.0062 


3 for Symanzik action at £r = 3. 




Symanzik Action 


at £r = 3 


P 


V 


^Nobub 


10.0 


1.0341 ± 0.0146 


1.0426 ± 0.0058 


8.0 


1.0260 ± 0.0150 


1.0361 ± 0.0042 


6.0 


1.0520 ± 0.0200 


1.0667 ± 0.0015 


4.5 


1.1482 ± 0.0317 


1.1331 ± 0.0064 



Table 4 

r\ for Iwasaki action at £r = 2. The data for (3 = 3.5 with * are calculated on the 
16 3 x 32 lattice to study the size dependences. 

Iwasaki Action at £r = 2 



p 


V 


oSub 


10.0 


0.9811 ± 0.0030 


0.9742 ± 0.0033 


6.0 


0.9831 ± 0.0037 


0.9784 ± 0.0032 


4.5 


0.9755 ± 0.0083 


0.9776 ± 0.0044 


4.0 


0.9806 ± 0.0074 


0.9782 ± 0.0039 


3.5 


0.9761 ± 0.0105 


0.9767 ± 0.0049 


3.05 


0.9911 ± 0.0182 


0.9881 ± 0.0060 


2.5 


0.9998 ± 0.0045 


0.9837 ± 0.0074 


3.5* 


0.9803 ± 0.0070 


0.9891 ± 0.0036 



in perturbative and numerical results. However, the slope of 77 becomes steeper 
for the numerical results. 
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Table 5 

77 for Iwasaki action at £r = 3. 



Iwasaki Action at £r = 3 



p 


V 


N oSub 


10.0 


0.9714 ± 0.0054 


0.9647 ± 0.0006 


6.0 


0.9569 ± 0.0041 


0.9554 ± 0.0026 


4.0 


0.9700 ± 0.0118 


0.9645 ± 0.0063 


3.5 


0.9715 ± 0.0160 


0.9708 ± 0.0031 


3.05 


0.9725 ± 0.0120 


0.9776 ± 0.0037 


2.56 


1.0011 ± 0.0138 


1.0011 ± 0.0071 



Table 6 

77 as a function of (3 for DBW2 action at £r = 2. 

DBW2 action at £ R = 2 



Table 7 

77 for DBW2 action 



p 


77 


oSub 


2.5 


0.9084 ± 0.0090 


0.8626 ± 0.0025 


1.6 


0.9011 ± 0.0082 


0.8616 ± 0.0018 


1.4 


0.8917 ± 0.0122 


0.8623 ± 0.0024 


1.2 


0.8882 ± 0.0115 


0.8673 ± 0.0032 


1.1 


0.8868 ± 0.0144 


0.8753 ± 0.0030 


1.0 


0.8781 ±0.01069 


0.8817 ± 0.0092 


at Cr 


= 3. 






DBW2 action at £ R = 3 


P 


77 




1.4 


0.8283 ± 0.0189 


0.8082 ± 0.0046 


1.2 


0.8157 ± 0.0252 


0.8070 ± 0.0055 


1.1 


0.8122 ± 0.0230 


0.8210 ± 0.0076 


1.0 


0.8123 ± 0.0235 


0.8262 ± 0.0101 



In this case the tadpole improved one-loop perturbation calculation (boosted 
perturbation) [22,23] reduces the discrepancy. It may be to replace (5 in Eq. 4 

by = PyJW a (l, l)W t (l,l): 

V&0) = l + % , ^ , = • (12) 
P ^.(1,1)^(1,1) 



In this formula, since W s (l, 1) and W t (l, 1) decrease as P decreases, the P 
dependence of 77 is more enhanced. The use of Eq.12 improves the agreement 
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17: Symanzik and Standard Action 



1.30- 




1.25- 


l 




J 


1.20- 


_L 


1.15- 


1 


1.10- 


f 


1.05- 





1.00- 


1 1 1 



• v 2 
■ v 5 



H 2 s (Pert) 



V 2 (b-pert) 
V 3 S (Pert) 



-°--y 3 (b-pert)) 



^2 

7j^ and (Pert) 




Fig. 4. 77s for Symanzik (rj Sz ) and standard (j] Stand ) actions. For standard action, 
we show the results obtained by Klassen[2]. (5 of standard action is shifted by Eq.5, 
in order to compare the 77 at the same lattice spacing a. The perturbative results 
are also shown for comparison with numerical results. The 77(6 — pert) represents 
the boosted perturbation result(Eq. 12). 



between perturbative and numerical results a little as shown in the Fig. 4. 

In this figure, we also show the result from the standard plaquette action[2]. 
In order to compare 77 at the same lattice spacing, we have shifted f3 of standard 
action to that of Symanzik action using the asymptotic scaling relation Eq. 
5. Because the two estimations coincide for Symanzik action for Pent at the 
N T = 8 lattice in Section 2.1, in these regions the asymptotic scaling relation 
Eq.5 may be satisfied at least approximately for these two actions. It is found 
that the qualitative behaviors are the same, although the slope becomes more 
gentle for Symanzik action. 



3.3 Iwasaki Action 



Results for Iwasaki action are shown in Fig. 5. 772 exhibits a shallow dip 
around /3 ~ 4.5 and then increases with decreasing f3. 773 shows similar behav- 
ior, but the position of the dip moves to (5 ~ 6.0. Both 772 and 773 stay close 
to unity in a wide range of j3 for (3 > 2.5. The deviation from unity is more 
enhanced for 773. 
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Fig. 5. rj for Iwasaki action. 

In the continuum limit, the 77 parameter approaches one. The improved 
action has an rj value that remains one in a wide range of /3; 2.5 < /3. The 
deviation from unity is within 4% when = 2, and 3, therefore until precision 
simulation of a few % accuracy is required, detailed calibration of 77 is not 
necessary. This is a good property for simulations. 

For Iwasaki action, the one-loop perturbative calculation predicts a mono- 
tonic decrease in rj as (3 decreases[4], as shown in Fig. 5. The numerical results 
are qualitatively different from those of the one-loop perturbation. 



3.4 DBW2 Action 



Results for DBW2 action are shown in Fig. 6. r] stays almost constant; 7] 2 ~ 
0.9 in the range 1.0 > (3 > 2.5 and r] 3 ~ 0.81 in the range 1.0 > (3 > 1.4. 
The deviation from unity is not small in this case, but the flatness of rj as a 
function of (3 is again a good property for numerical simulations. The rough 
calibration of 77 at a few f3 points is sufficient to obtain a reasonable estimation 
of £b. As in the case of Iwasaki action, the numerical results show qualitatively 
different behavior from the one- loop perturbative ones [4] (see Fig. 6). 
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11 

1.0-. 




Fig. 6. rj for DBW2 action. 
4 Discussion and Conclusions 

In this work, we studied the global structure of rj as a function of /3 and C\ 
for the gauge action given in Eq.l. Overall effects of the improved actions on 
7] are summarized as follows. The plaquette term in the action makes the r\ 
parameter increase monotonically as (3 decreases, while the rectangular terms 
with C\ < make 77 decrease. 

At C\ = —1/8, the effects of rectangular loops are not so large, and the 
slope of i] is smaller than that in the case of standard action. As a result, at 
the same lattice spacing, the effects of the quantum correction are reduced in 
Symanzik action. In Symanzik and standard actions, the 77 dependence on /3 is 
qualitatively the same between perturbative and numerical results, but slopes 
are steeper for numerical results. 

At —C\ = 0.331, the increase in rj due to the plaquette term and its sup- 
pression by the rectangular loops are almost in balance in a wide range of /3, 
2.5 < j3, and 77 stays close to one. However, the detailed behavior of these 
effects depends on (5 and £r. r] 2 and 773 decrease from (5 = 2.5, and exhibit 
shallow dips at f3 ~ 4.5 and f3 ~ 6.0, respectively. These behaviors are quali- 
tatively different from the prediction of one-loop perturbative calculation. 

At —C\ = 1.4088, the effects of the rectangular loop become stronger than 
those of the plaquette term, and 77 becomes less than one. It is almost inde- 
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pendent of (3 in the range 1.0 < (3 < 2.5. This is again qualitatively different 
from the perturbative result. 

As increases, the effects of the two terms in the action are more enhanced, 
and the balance becomes more subtle. 



In the continuum limit, r] approaches unity. Then Iwasaki action is close 
to the continuum limit in the region 2.5 < (3 for £r = 2 and 3. Particularly 
around (3 ~ 2.5, r\ is close to unity. This means that the calibration of i] is not 
necessary until high precision simulation is carried out. 

In the case of DBW2 action, r\ is not close to one. Then, as far as i] is 
concerned, it is not close to the continuum limit in this (3 region. However, 
i] is almost independent of (3. This is a good property for the simulation of 
physical quantities on anisotropic lattices, because the rough calibration of r\ 
is sufficient to determine in this action. 

For Symanzik action, the deviation of i] from unity is remedied compared 
with standard action, but the slope of i] is steep, and it becomes ~ 10% at 
around (3 ~ 5.0. Therefore detailed calibration of r\ is necessary. 

For the (3 and £ ranges that we have studied, the differences between rj 
and rj Nosub are small for all the improved actions. For Symanzik and Iwasaki 
actions, the difference is ~ 1%, and for DBW2 action, it is a few %. Therefore 
it is safe to use rj Nosub except in the case of very precise simulations. This is 
good news, because the calculation of 77 requires much CPU time. 

Further data on r] for larger £r and smaller (3 will be reported in the forth- 
coming publications, because the calculation of r] at smaller (3 and larger £r 
requires much more CPU time. 
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A Optimal Radius of Integration For Iwasaki and Symanzik Ac- 
tions 



If R is an external source field for link variable U, the link integration of U 
is given by 

_ 1 dZ(R) _ J D[U]Uexp(TrjRW + Utf)) 
[ ' Z dW J D[U]exp(Tr(RW + UR))) 1 } 



where Z(R) is expressed by the modified Bessel function Ii [9,10]. 



and 



P(x) 



X 



Q = 2Re(det(R)), 

P(x) = l + xTr(RR^) 

+\x 2 [(Tr(Rrt)) 2 -Tr((R&) 2 



x 3 det(RR*). 



(A.3) 



Similarly dZ[R)jdR is written by the modified Bessel function I x and I 2 [9,10]. 

dZ(R) r dx , rCl \ r . . dQ r dx e xQ r , n ,dP(x) , . , x 
— i_Z = i xe xQ -I i(2z)-^r + * — I 2 (2z) pr-V (A.4) 



The path of the integration is a closed circle on the complex plane x. In 
principle it is arbitrary, but numerical integration requires adequate radius. In 
the case of the standard action, the adequate radius was studied by Scheideler 
[24]. 

The arguments of the modified Bessel functions become rather large and we 
apply asymptotic expansion. In this article we use Simpson method for the 
numerical integration, and search for the region of r and the number of the 
division N where (U) is stable against the change of r. 

An example of the r dependence of a (U) is shown in Fig.A.l. It is found 
that when the number of the devision is N = 100, there appears some spurious 
plateaus, which disappear when N = 400. However there is a region of r where 
(U) is stable in the change of N, which is the optimal region of integration for 

= 100. The optimal region increases a little when A^ = 400. In this article 
we choose A^ = 100 and proceed to determine the optimal region of r (r opt ). 
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0.01 



Fig. A.l. Radius dependence of (U4(2, 2)). The integration is carried out on a config- 
uration of Iwasaki action, fully thermalized at (3 = 3.05 and £b = 2.0. An integrated 
link is located at the center of the lattice, directing in the 4-th direction. We have 
shown the (2,2) elements of the unitary matrix U. 



These plateaus shown in Fig. A.l are observed when Taylor expansions of 
the modified Bessel functions are applied. Then they are due to the difficulty 
in the numerical integrations given by Eq.A.2 and Eq.A.4. Therefore it is 
important to find the optimal r region. 

For some set of f3 and we have obtained the minimum and maximum 
of r opt as shown in the Fig. A.l, for space-like and time-like links separately. 



When £ > 1.0, r opt of time-like links 



r opt\ 



is smaller than that of space-like 



links (r° pt ). Examples of the difference is shown in the Fig. A. 2. It seems that 
the difference becomes larger as (3 and increase. The similar properties are 
observed in the case of standard action. 



We proceed to the parameterization of the r op '(/3,£). The (3 and £ ranges 
are 2.0 < (3 < 6.0, 1.8 < £ B < 6.5 for Iwasaki action and 4.5 < (3 < 8.0, 1.7 < 
£b < 5.8 for Symanzik action. The r° pt ace (min) and r° pt ace (max) are shown in 
the Fig. A. 3. They decrease with (3 and £ and seems to be parametrized as 



r opt = a * exp(-b(3 - c£ B ) 



(A.5) 
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r 

0.1 -- 



0.01 



1E-3 



2 3 4 a 5 6 7 



Fig. A. 2. Examples of the difference between r°£* ce and r^j e for Iwasaki action at 
Cn = 2.0. 




Fig. A. 3. A compilation of R m in, (a) and i? max (b) of Iwasaki action in the range 
2.0 < < 6.0 and 2.0 < £ B < 6.0. 

Then we define r opt (mid) as, 

\og(r opt (mid)) = (\og(r opt (min)) + \og(r opt (max)))/2. (A.6) 

and then fit them by Eq. A. 5. For log(r op '), it becomes a multiple regression 
fitting. The coefficients are determined by the least square method. The re- 
sults for Symanzik and Iwasaki actions are summarized in the table A.l. We 
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Dpt Iwasaki: £ =2.0 



t-direction 

V r 0F >ax) 

A r 0pt t (min) 
s-direction 



r 0pl s (max) 
r 0pt (min) 



W 

A 



Table A.l 

The multiple regression fit of r opt (mid). The 60 data and 29 data points are used 
to determine the coefficients a, b, c for Iwasaki and Symanzik actions respectively. 



Action 




a 


b 


c 


Symanzik 
(26data) 


r° pt (mid) 
r^iynid) 


0.5563 
0.06244 


0.5479 
0.4213 


0.5336 
0.6568 


Iwasaki 
(61data) 


r° pt (mid) 


0.08663 
0.01682 


0.5507 
0.5139 


0.4315 
0.5261 



have checked that r opt (mid) with the parameter given by table A.l is located 
between r opt (max) and r opt (min); namely it stays within the optimal radius 
of integration through out the data points. 

r opt region may change with the position of link on a configuration and also 
with configurations. The results shown in Fig.A.3 are obtained for a link at 
the center of the configuration in space and time directions, which are fully 
thermalized. However the fluctuation of the r opt region is not large. If we 
choose, r opt (mid), it has been in an optimal region of r for all link variables 
and configurations. 

Let us proceed to discuss the effects of the link integration method. In 
the case of improved actions, the number of link U which are simultaneously 
integrated in a Wilson loop becomes much smaller than the case of standard 
action, because in the case of improved action the background fields R of 
Eq. A. 2 extend wider range due to the 6- link rectangular loop in the action. 
Therefore the effect of link integration method is reduced in these cases, and 
is not effective for the calculation of smaller Wilson loops. 

The example of the suppression of the fluctuation is shown in the Fig. 
A. 4. The suppression of the fluctuation is impressive for W(6, 6) but not 
for W(A, 4). The similar properties are observed for the Symanzik action of 
W(8,8) and W(4,4) at rj = 4.5. The link integration needs much CPU-time, 
and therefore the link integration method is effective for the calculation of 
large Wilson loops in the confined phase or very close to the transition point. 
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